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ABSTRACT 

We report a Chandra ACIS-I observation of the dense globular cluster Terzan 5. The previously known transient 
low-mass x-ray binary (LMXB) EXO 1745-248 in the cluster entered a rare high state during our August 2000 
observation, complicating the analysis. Nevertheless nine additional sources clearly associated with the cluster 
are alsodetected, ranging from Lx(0. 5-2. 5keV) = 5. 6 x 10^^ down to 8.6 x 10^' ergs s"'. Their X-ray colors and 
luminosities, and spectral fitting, indicate that five of them are probably cataclysmic variables, and four are likely 
quiescent LMXBs containing neutron stars. We estimate the total number of sources between Lx(0.5-2.5keV) = 
lO"*^ and lO"'-' ergs s"' as 11. 4;!^ ^ by the use of artificial point source tests, and note that the numbers of X- 
ray sources are similar to those detected in NGC 6440. The improved X-ray position allowed us to identify a 
plausible infrared counterpart to EXO 1745-248 on our 1998 Hubble Space Telescope NICMOS images. This 
blue star (Fl 10W=18.48, F187W=17.30) lies within 0.2" of the boresighted LMXB position. Simultaneous Rossi 
X-ray Timing Explorer {RXTE ) spectra, combined with the Chandra spectrum, indicate that EXO 1745-248 is an 
ultracompact binary system, and show a strong broad 6.55 keV iron line and an 8 keV smeared reflection edge. 



Subject headings: X-rays: individual (EXO 1745-248) — 
globular clusters: individual (Terzan 5) - 

1. INTRODUCTION 

The high resolution of the Chandra X-ray Observatory has 
enabled astronomers to study the low-luminosity X-ray source 
populations in globular clusters in great detail. Combined X- 
ray, radio, and optical Hubble Space Telescope (HST ) stud- 
ies of the globular cluster 47 Tucanae have revealed quies- 
cent low-mass X-ray binaries (qLMXBs) that have not expe- 
rienced X-ray outbursts in the history of X-ray astronomy, cat- 
aclysmic variables (CVs) as X-ray luminous as any known in 
the field, flaring behavior from coronally active stellar binary 
systems, and predominantly thermal X-ray emission from mil- 
lisecond pulsai-s (MSPs) (Grindlay et al. 2001a, 2002). Similar 
populations have been uncovered in the globular clusters NGC 
6397 (Grindlay et al. 2001b), NGC 6752 (Pooley et al. 2002a), 
and u! Cen (Rutledge et al. 2001, Cool, Haggard & Carlin 
2002), while the luminosities and broad X-ray spectral types 
of these sources have allowed classification of sources in the 
more obscured cluster NGC 6440 (Pooley et al. 2002b) and 
M28 (Becker et al. 2003). 

The globular cluster Terzan 5 contains a known transient 
LMXB, EXO 1745-248, which was first detected in a burst- 
ing state by Hakucho (Makishima et al. 1981), and has been 
irregularly active since then (Johnston et al. 1995 and refs. 
therein). EXO 1745-248 is one of the few luminous globular 
cluster LMXBs not analyzed by Sidoli et al. (2001; hereafter 
SPOOl) or Parmar et al. (2001), who identify spectral distinc- 
tions between normal and ultracompact LMXBs. Ultracompact 
LMXBs, defined as having periods less than 1 hour, are thought 
to possess a degenerate helium white dwarf secondary. Deutsch 
et al. (2000) remark upon the overabundance of ultracompact 
LMXBs in globular clusters, and speculate that the short peri- 
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ods may be due to dynamical effects in globular clusters. 

Terzan 5 also contains two identified MSPs, with many addi- 
tional MSPs probably making up the extended steep-spectrum 
radio source at the cluster core (Lyne et al. 2000; Fruchter & 
Goss 2000). Terzan 5's high central density and large mass 
make it a rich target for studies of binary systems, but its high 
reddening and severe crowding make optical and even infrared 
observations extremely difficult. The deepest infrared survey 
of Terzan 5 was performed with the HST NICMOS camera by 
Cohn et al. (2002; hereafter CLG02). The extreme reddening 
indicates that the cluster parameters are best determined in the 
infrared. NICMOS observations produced the first CMDs of 
Terzan 5 to reach the main-sequence turnoff (CLG02; see also 
Ortolani et al. 2001). Therefore we utilize the new cluster pa- 
rameters derived by CLG02 in our analysis, particularly the red- 
dening, distance, core radius, and radial star-count profiles. The 
attempt by Edmonds et al. (2001, hereafter EGCOl) to iden- 
tify EXO 1745-248 and the eclipsing MSP through time-series 
variability analysis and color information did not uncover any 
promising candidates due to the crowding and overlapping Airy 
profiles, although it did identify two truly variable stars, one of 
which was shown to be an RR Lyrae variable. 

This paper is organized as follows. Section 2.1 describes the 
observations we used. Section 2.2 explains our methods for de- 
tecting sources and performing an astrometric correction based 
upon identification of serendipitous sources. Section 2.3 de- 
scribes our search for an infrared counterpart for EXO 1745- 
248, extending the work of EGCOl using our Chandra posi- 
tion. Section 2.4 classifies the faint X-ray sources, and quanti- 
fies our detection incompleteness (due to the outburst of EXO 
1745-248). Section 2.5 attempts simple spectral fits to the faint 
sources, while section 2.6 examines the simultaneous RXTE 
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and Chandra spectra of EXO 1745-248. Section 3.1 compares 
our spectral analysis of EXO 1745-248 to other observations of 
LMXBs. Section 3.2 recalculates the central density and col- 
lision frequency of Terzan 5, while section 4 summarizes our 
conclusions. 

2. ANALYSIS 
2.1. Description of Obse rvations 

The Chandra X-ray Observatory observed the globular clus- 
ter Terzan 5 on July 24, 2000, for 45 ksec (05:46 to 18:22 TT), 
and on July 29, 2000 for 5 ksec (00:56 to 09:20 TT), with the 
AClS-1 instrument at the focus. Due to an error in the obser- 
vation upload, both exposures were performed in 1/8 subar- 
ray mode (as intended for the shorter observation only), with 
frame times of 0.841 s (longer exposure with more chips) and 
0.541 s. EXO 1745-248 entered an outburst during July 2000, 
its 2-10 keV flux varying between 54 and 600 mCrab during 
July and August, approaching its Eddington limit at maximum 
flux (Markwardt et al. 2000a, 2000b). Rossi X-ray Timing Ex- 
plorer (RXTE ) All-Sky Monitor (ASM) observations (results 
provided by flie ASM/RXTE team^) show that the 2-10 keV 
countrate from EXO 1745-248 on July 24, 2000, was 5.3±0.9 
cts s"' (~ 72 mCrab) during our longer Chandra observation, 
rising to 9.8±0.9 cts s~^ in the second observation (Fig. 1). 

Unfortunately, the Chandra observation of EXO 1745-248 
was not optimized to study such a bright object (Fig. 2). The 
intense photon flux led to severe pileup, which occurs when 
two or more photons landing in the same or adjacent pixels be- 
tween frames are recorded as a single event. Pileup can in- 
crease the energy of recorded events, or by altering the grade 
of the recorded event to a "bad" grade, cause the rejection of 
the event either before or after telemetry to the ground (See the 
Chandra Proposer's Observatory Guide, v.5, chapter 6). In this 
case the pileup was severe enough to cause the pre-telemetry re- 
jection of all events recorded within ^ 1 " of the LMXB position. 
The LMXB's X-ray halo, by increasing the local background, 
greatly degraded our sensitivity to faint cluster sources. (The 
halo is due to a combination of dust grain scattering and the in- 
trinsic breadth of the Chandra mirrors' point spread function.) 
Nevertheless, the spectacular resolution of Chandra did allow 
us to identify additional point sources up to 10^ times fainter 
within 10" of EXO 1745-248. The readout streak (caused by 
photons from EXO 1745-248 arriving during the readout of the 
CCD) also degraded our survey, but we were able to extract 
a useful spectrum of EXO 1745-248 from the readout streak. 
We also analyzed a simultaneous Rossi X-ray Timing Explorer 
(RXTE ) pointed observation from the HEASARC archive^ (on 
July 24, 2000, 15:15 to 16:16 TT), for broad spectral coverage 
of the outburst of EXO 1745-248. The HST NICMOS data we 
used to search for a possible infrared counterpart to EXO 1745- 
248 are described in EGCOl and CLG02. 

2.2. Detection and astrometry 

We used the ClAO software package-' to search for point 
sources, produce hardness ratios, and extract spectra and 
lightcurves. We reprocessed the two observations to remove the 

' See http://xte.mit.edu/. 

^ Available at http://heasarc.gsfc.nasa.gov/docs/corp/data.html. 

^ Available at http://asc.harvard.edu/ciao/. 

GSC 2.2 absolute astrometric errors, when compared to the international 
gsss.stsci.edu/gsc/gsc2/GSC2home.htm. 



pixel randomization added in standard processing, and merged 
the two observations. No periods of high background flaring 
were observed. EXO 1745-248's X-ray halo displays an x-ray 
color (defined here as 2.5 log [0.5-1.5 keV counts/1.5-6 keV 
counts], following Grindlay et al. (2001a)) of -2.9, harder than 
most known faint globular cluster sources. Therefore, we se- 
lected a soft band, 0.5-2 keV, to search for point sources. We 
selected a 2.5 x lO'' pixel (1.7 arcmin^) square region including 
the cluster and ran WAVDETECT with the significance thresh- 
old set to give false positives at the rate of 10"'*. WAVDETECT 
indeed found two spurious sources (identified as spurious by 
eye, and by not appearing in more than one energy band) far 
from the cluster, as expected. Within three optical core radii 
(24") we identify nine point sources besides EXO 1745-248, 
each confirmed by visual inspection (see Figure 2) and with 
significance above 2.8 cr. We name the nine additional real 
sources in Table 1 with both positional names and simple ref- 
erence names W2-W10 (used in the rest of this paper). Outside 
the globular cluster, we binned the remaining data into 1 arc- 
second pixels and searched for back- or foreground serendip- 
itous sources (setting the significance threshold to a false rate 
of 10"*). We identified four sources that also pass visual in- 
spection on the active portions (36 arcminutes^) of the AClS-1 
array. As they are all more than 2' (2.5 half-mass radii) from 
the cluster center, we hereafter identify them as "serendipitous 
sources". The ROSAT source S2 from Johnston et al. (1995) 
was not included in the field of view. Inspection of the locations 
of the known millisecond pulsars (A and C; Lyne et al. 2000) 
show no evidence of X-ray emission (using the astrometry be- 
low). The upper limits are 3 and 5 counts respectively in the 
0.5-2 keV band, giving < 2 x 10^' and < 3 x 10^' ergs s"' 
for a power law of photon index 2. These are well above the X- 
ray luminosities of all identified millisecond pulsars in 47 Tuc 
(Grindlay et al. 2002), NGC 6397 (Grindlay et al. 2001b), and 
NGC 6752 (D'Amico et al. 2002), but below flie luminosity 
{Lx = 1.1 X 10^3 ergs s"') of the MSP B1821-24 in M28 (Dan- 
ner et al. 1997). Soft thermal spectra give even less constraining 
limits, due to the high extinction. The two variables identified 
by EGCOl are also not seen, as expected based on their RR 
Lyrae and probable eclipsing blue straggler identifications. 

Inspection of a Digital Sky Survey image reveals that one 
of the serendipitous sources has a probable optical counter- 
part. According to the Guide Star Catalog 2.2, a star with 
V=13.56, R=13.22 is located only 0"36 from flie position de- 
rived by WAVDETECT for the X-ray source that we name 
CXOU J174803. 3-244854. This is consistent with the 0"6 ab- 
solute astrometry (90% conf . radius) reported by Aldcroft et al. 
(2000) for Chandra . We estimate a probability of a positional 
coincidence within 0"5 of one of the 15 brightest stars on the 
7' by 7' survey plate with any of the four serendipitous X-ray 
sources of 8 x 10"^ Another X-ray source (CXOU J174812.6- 
2448 1 1 . 1) has a faint star (R=16.8) Of 14 away (when the frame 
is shifted to match star 1). We estimate the chances of a star this 
bright or brighter landing within 0."5 of one of four Chandra 
sources as 4 x 10""'. The other two X-ray sources show no stars 
in the GSC 2.2 Catalog within 3". We shift our Chandra frame 
by -i-0?004 in RA, and H-0."396 in Dec, to match flie weighted 



celestial reference frame, are of order 0.3-0.75 arcseconds; see http://www- 
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GSC 2.2 star positions."* 

The V and R magnitudes of star 1 suggest an F star, and 
the X-ray (0.5-2.5 keV) to V-band flux ratio (V-band flux de- 
fined as IQ-O'^v-S '^'i ergs cm^^ s'^) is 2 x IQ-^, on the high end 
of values for nearby F stars (Hiinsch et al. 1999). Star 2 has 
no color information in the GSC 2.2 Catalog. Assuming V- 
R'~1.5 (appropriate for M2 stars) gives a flux ratio of 1 x 10"^, 
which is common among M stars (Hiinsch et al. 1999). The 
other two sources may be background AGN, or (perhaps more 
likely) CVs in the galactic bulge. Both show an Xcolor (2.5 
log [0.5-1.5 keV cts/1.5-6 keV cts]) of -1.5, indicating strong 
absorption. Simple spectral fitting of the brighter source gives 
Nh = 2.8^i 2 X 10^^ cm"^ for a powerlaw of photon index 1.7, 
as typical for AGN. A 10 keV bremsstrahlung spectrum (as for 
a bulge CV) would have A^^ = l.Tf^ x 10^^ ^.^-2^ j^^^ 2-10 
keV log N-log S relation of Giacconi et al. (2001) suggests 
that 1-3 AGN may be expected at these observed flux levels 
in our ACIS-I subarray field (although the location of detected 
sources near to the aimpoint suggests that our sensitivity is not 
uniform across the field). Recent results from galactic bulge 
surveys (Grindlay et al. 2003) indicate that bulge CVs or com- 
pact binaries outnumber AGN at these flux levels, suggesting 
that these two sources may be bulge CVs. 

We found the position of EXO 1745-248 by centering the 
symmetric "hole" caused by extreme pileup and the readout 
streak. We estimate our error in this determination at 2/5 of 
a pixel, or 0."2. We present the positions (in the GSC 2.2 
frame) of EXO 1745-248, the nine additional globular cluster 
sources, and the four serendipitous sources (with relative po- 
sitional errors from WAVDETECT) in Table 1, along with the 
background- subtracted counts in three bands. 

2.3. Identification of plausible infrared counterpart to EXO 

1745-248 

Our refined position for EXO 1745-248 aUowed us to under- 
take a careful search of the small Chandra error circle in the 
June 1998 Hubble Space Telescope NIC2 FllOW and F187W 
data, for which the photometric analysis is described in EGCOl. 
We estimate an error of 0."1 in the shift to the GSC 2.2 frame, 
based on the scatter of the two stars used in the boresighting. An 
additional error of 0"37 (following EGCOl) is incurred in the 
shift between the GSC 2.2 and HST NICMOS frames. Com- 
bining these with the 0"2 centering positional uncertainty gives 
an error for the position in the NICMOS frame of 0."4. This er- 
ror circle is fortuitously free of bright red giants. Error circles 
for six other X-ray sources in the cluster lie wholly or partially 
within the NIC2 field of view. 

Comparison of the FllOW and F187W images with a ratio 
image (Figure 3) shows a significantly blue faint star (hereafter 
star A) within the Chandra error circle, only 0."2 from the cor- 
rected position of EXO 1745-248. A combination of this star 
and another star just above it (Star B) was shown in the cleaned 
CMD of EGCOl (their figure 8) as one of only two stars that 
lie more than 3c7 to the blue from the distribution of stars in 
that color-magnitude diagram. Our attempts at separate PSF- 
fitting using ALLSTAR in DAOPHOT encountered two prob- 
lems: (1) we were unable to fit stars A and B simultaneously 
and (2) the sky-fitting component in PHOT and ALLSTAR gave 
local sky (background) values that were too high, because stars 
A and B are surrounded by a ring of bright stars. The latter is a 
formidable problem in this crowded field, especiaUy in F187W 
where the background contribution of neighboring red giants is 



greater than in Fl lOW. 

We therefore adopted a customized procedure to estimate the 
FllOW and F187W magnitudes for stars A and B. For sky de- 
termination we adopted the mean of the two 2-by-2 pixel re- 
gions having the lowest intensity values within an 8 pixel radius 
of star B. Using the FllOW and F187W images, and the ratio 
image, we determined the positions of stars A and B. Using 
PSF models and the DAOPHOT program ADDSTAR to create 
models for stars A and B we subtracted these models from the 
two images to produce residual images, while keeping the po- 
sitions of stars A and B fixed. Since no other close-by stars 
are obviously present, stars A and B are accurately modeled 
if the mean aperture counts (after subtracting the sky) at their 
positions in the residual image are equal to zero. We iteratively 
varied the magnitudes of these two stars until this goal was met. 
The results of these fits are shown in the lower right-hand cor- 
ner of Figure 3. Clearly the subtractions of stars A and B were 
successful, given the relatively smooth residuals. 

We performed aperture photometry on the residual images to 
produce photometry that is consistent with Figure 8 of EGCOl . 
Photometry for stars A and B were measured by subtracting 
off stars B and A respectively, and the results are shown in 
Figure 4. Variables 1 and 2 from EGCOl are labeled with 
crosses, and stars within 1" of X-ray sources are circled. The 
main-sequence tumoff (not reached at the cluster core) is at 
F110W« 20 (CLG02). Note the apparently blue color of star 
A (~ 2.5cr from the subgiant sequence) lying near the top of 
the blue straggler sequence. The error in the color of star A 
is dominated by the error in the F187W magnitude. The high 
sky background in F187W caused the original F187W magni- 
tude determination for the combination of stars A and B to be 
too low, and thus shifted the combination bluewards. Our cus- 
tomized procedure thus gives redder colors for both stars A and 
B than the original determination for the combination. 

Given the blue appearance of star A in both the normal and 
residual images, the reasonable evidence for blue colors in our 
photometric analysis, and its proximity to the center of the X- 
ray error circle, star A is a reasonable candidate for the IR coun- 
terpart to EXO 1745-248. No evidence is seen for statistically 
significant variability of star A, but the faintness of star A and 
its strong blending with star B limits the depth of variability 
searches. Only the FllOW images are useful for time series 
analysis, and only sparse coverage was obtained here over 1.16 
days. Images in the V and / passbands with the Advanced Cam- 
era (in High Resolution mode) on HST should do a much better 
job of separating stars A and B and would provide much bet- 
ter constraints on variability if taken at several different epochs 
near outburst and quiescence. 

Since the F187W band does not correspond to any ground- 
based filter, standard calibrations for it do not exist. We ap- 
ply rough transformations to estimate interesting quantities. 
We estimate the unreddened J magnitude from the Fl lOW and 
F187W magnitudes by the equation 

7o = /^110W + 0.033-0.178(fllOW'-F187W')+A., 
derived assuming that the spectral slope is constant over this 
range. Using Ay/Ay=0.282 and Av=6.15 (CLG02), we find 
70=16.40, My=1.7. We calculate the color (F110W-F187W)o= 
-0.17 ± 0.3. Assuming that the photometry is accurate, three 
alternatives to an LMXB nature remain for this star: it could 
be a blue horizontal branch star, it could be a blue straggler, or 
it could be a foreground star. We plot simple interpolated blue 
straggler and horizontal branch sequences in Figure 4. Star A's 
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position may be consistent with that of a blue horizontal branch 
star, although a blue horizontal branch has only been hinted at 
in Terzan 5 (CLG02). It is also consistent with an extremely 
heavy (more than twice the tumoff mass) blue straggler, or a 
blue straggler binary (as variable 2 from EGCOl, plotted as an 
"X" next to star A in figure 4, may be). The other extremely 
blue star seen in Figure 4 is probably a foreground star, which 
also cannot be excluded as a possibility for star A. However, 
the chance of one of the five bluest stars in the 19'.'2 square 
NIC2 FOV landing in the 0."4 error circle by chance is only 
7 X 10"^. None of the stars in the other X-ray error circles (cir- 
cled in Figure 4) display any unusual colors or variabiUty. This 
is expected based upon the magnitudes of CVs and quiescent 
LMXBs, in the field and other globular clusters, which very 
rarely reach the turnoff magnitude (Grindlay et al. 2001a, b, 
Pooley et al. 2002a). 

We compare this suggested counterpart to other low-mass X- 
ray binaries observed in the infrared by calculating the ratio of 
X-ray to infrared luminosity. Using the ASM lightcurve in fig- 
ure 1, we see that the average ASM countrate for the 1 1 days 
leading up to the HST observation was 1.5 ± 0.2 cts s~', com- 
pared to 5.3±0.9 cts s"^ at the time of the Chandra observa- 
tion. The high source confusion near the galactic center for 
the ASM makes this data relatively noisy, so we take 1.5 cts 
s~' as an upper limit on the X-ray flux. Assuming the same 
spectrum as we find in 2000 (see Sect. 2.5), we derive an un- 
absorbed 0.5-10 keV X-ray flux of 6.1 x lO"'" ergs cm'^ s'^ 
and an unabsorbed J-band flux of 1.2 x 10"'^ ergs cm"- s"', 
thus log(Lx/Lj)<3.7. We note that van Paradijs and McClin- 
tock (1995) find log(Lx(2-ll keV)/Lo/>/ (3000-7000 A))~ 2.7 
for typical bright LMXBs. 

Rough (./ - K\) colors have been reported for the soft X-ray 
transients in outburst QX Nor, or X1608-52 (0.5 < {J-K)o < 
-0.4, Wachter 1997), and SAX J1808.4-3658 (0.1 ±0.1, Wang 
et al. 2001) and flie faint LMXB XI 323-619 (-0. 1 ± 1 .4, Smale 
1995). Our value of (Fl 10W-F187W)o - -0.17 ±0.3 for star 
A, implying (J - K)o ~ -0.3, is consistent with this broad 
range. Infrared colors observed by Wachter (1998) give an even 
broader range (-0.4 < (J-K)o < 1.41), as some LMXBs are 
dominated by the disk light in the IR, while others are domi- 
nated by the light of the secondary. This makes identification 
of LMXBs by infrared colors very difficult, and leaves our iden- 
tification uncertain. 

2.4. X-ray Classification and Completeness 

We extracted the numbers of counts within 1" radius of the 
sources (containing 90% of source flux at 1.5 keV^) in the 0.5- 
1.5, 0.5-4.5, and 1.5-6 keV bands. We then subtract 1/3 (ratio 
of areas) of the counts from background l"-2" annuli around 
these sources for each band (except where the annuli would 
overlap with other sources, where we offset the extraction re- 
gion). The results are listed in Table 1. Constructing an X-ray 
CMD (as in GHEOla; Xcolor is defined as 2.5 log[0.5-1.5 keV 
counts/1.5-6 keV counts], while the vertical axis is log[0.5-4.5 
keV counts]), we are able to identify two populations of faint 
sources with distinct X-ray colors. Using an Nh of 1.2 x 10^^ 
cm"^ (derived from CLG02's estimate of Av=6.75, using Pre- 
dehl & Schmitt's (1995) conversion factor), we apply a uniform 
shift to tiie coordinates of the CMD as Pooley et al. (2002b) do 
for the moderately reddened globular cluster NGC 6440. The 

^ CXC Proposer's Observatory Guide, http://cxc.harvard.edu/udocs/docs/docs.html 
* Available at http://space.mit.edu/ASC/MARX. 



identified dim sources (and the color of the LMXB halo, which 
determines the background) are plotted in the upper panel of 
Figure 5. The upper and right axes show the raw (absorbed) 
Xcolor and counts, while the lower and left axes represent the 
corrected (unabsorbed) Xcolor and counts. (We note that the in- 
terstellar medium decreases soft fluxes more than hard fluxes. 
Thus this simple shift of axes should be used with caution.) 
Another version of the same diagram, but plotting luminosities 
derived from spectral fits (section 2.5 below) against corrected 
Xcolor as above, is shown in the lower panel of figure 5. Errors 
on the vertical axis are derived from errors in the counts, and do 
not include spectral uncertainties. This version of the diagram 
(inspired by Becker et al. 2003), is more useful for comparing 
different clusters, but is dependent on spectral and reddening 
assumptions. The luminosities of the two populations roughly 
agree with the range observed for their respective counterparts 
in 47 Tuc and NGC 6440. This suggests that the soft population 
(W2, W3, W4, W8) should be identified with qLMXBs contain- 
ing neutron stars and that the harder population is composed of 
bright CVs. 

We estimate our incompleteness, due to the X-ray halo 
and streak from flie LMXB EXO 1745-248, by flie use of 
artificial-star tests. We generated uniform grids of artificial 
point sources at different luminosities using the MARX soft- 
ware^ and merged the simulated data with our 0.5-2.0 keV mas- 
ter image to see if we could detect the fake sources. We scaled 
the average (Poisson-distributed) counts of the artificial sources 
to the average detected counts of real sources for which XSPEC 
fitting gave the target luminosity (see section 2.5 below). E.g., 
for Lx = 10^^ erg sources, the average number of counts per 
simulated source was 20, and for Lx = 3 x 10^^ erg s"' it was 
45 counts. (The counts-luminosity relation varied depending 
on the kind of sources prevalent at different luminosities.) We 
performed the tests for artificial source 0.5-2.5 keV luminosi- 
ties of 10^2, 1.7 X 10^2, 3 X 10^2, and 10^^ ergs s"\ assuming 
spectra similar to the average of our detected sources. We sep- 
arated the artificial sources by 5" so that each source would 
suffer only the crowding that exists in the real data. We merged 
grids of point sources with the Chandra image four times for the 
10^^ and 1 .7 X 10^^ tests, and three times for the other two tests, 
shifting the grids by several arcseconds between each trial. We 
then apphed WAVDETECT to each image, keeping the same 
parameters as in section 2.2. We counted the number of artifi- 
cial sources that landed within one optical core radius, between 
one and two core radii, and between two and three core radii for 
each trial, and the number of sources detected in each region, 
for a total of ~1000 point-source trials. We note that this ap- 
proach takes crowding into account, in that the test sources will 
be crowded by any real undetected sources. Parts of the 1 -2 core 
radii annulus and 2-3 core radii annulus were imaged only in the 
5 ksec exposure and not in the 42 ksec exposure. The 2-3rc bin 
also includes a small area that is completely unexposed (due to 
the use of the subarray). We scaled the MARX point source 
luminosities appropriately and performed separate trials for the 
low-exposure regions (where sources below 3 x 10^^ ergs s~' 
were never detected), and then included these results in the to- 
tal completeness calculations. This aUowed an incompleteness 
factor to be determined for each radial annulus for each of four 
luminosity bins, which we created by equally dividing (in a log- 
arithmic sense) the ranges between our test luminosities. 
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The number of real sources detected was too small to ef- 
fectively determine both a radial distribution and a luminosity 
function (1 source within 1 5 within 1-2 r^, 3 within 2-3 
Tc). Therefore we assumed a luminosity function of the form 
dN oc L~''dlnL, L > L^in following Johnston & Verbunt (1996). 
We fixed 7 at 0.5 first, as found for sources between 10^^ and 
10^^ ergs s-i in NGC 6440 (Pooley et al. 2002b) and the gen- 
eral globular cluster system (Johnston & Verbunt 1996), and 
then allowed 7 to vary within the range 0.29 to 0.78 found by 
Pooley et al. (2002b). The uncertainty in our incompleteness 
due to different possible values of 7 was small compared with 
errors due to small-number statistics. For each annulus, we de- 
termined an average incompleteness factor using the results of 
our testing and the weights derived from the luminosity func- 
tion for the luminosity bins as specified above. 

To calculate how many sources were missed, we take the 
number of sources detected in each annulus and apply the la 
errors for small numbers given by Gehrels (1986). These num- 
bers and errors are multiplied by 1/f, where f is the complete- 
ness fraction, to obtain the average density of sources in each 
annulus. Then we calculated the expected numbers of missed 
sources, using the estimated density and again applying small- 
number errors. We ignore possible systematics in the method, 
since the random errors are so large. Thus, above 10^^ ergs s~\ 
we find that within 1 rc we predict 0.8^q| additional sources, 
within 1-2 we predict 1-5^^ 3 additional sources, and within 
2-3 rc we predict 1 . lilij Q sources, giving a total probable num- 
ber of 11 .4;!^ g sources between 10^^ and 10^^ ergs in Terzan 
5. We note that this source population is very similar to the 1 1 
sources in this flux range identified in NGC 6440 (Pooley et 
al. 2002b), as expected considering the similar cluster parame- 
ters (as derived in section 3.2 below). 

The radial distribution we infer can be compared with that 
predicted by a "generalized King model" (Lugger, Cohn, & 
Grindlay 1996) where the surface density of sources varies as 

where ax = -Oq- 1) and q = the ratio of the average 

X-ray system mass to the mass of the stars that dominate the 
density profile (i.e., tumoff-mass stars for King model clusters). 
The soft sources in 47 Tuc (thought to be largely MSPs) have 
ax = -3.6 and the CVs in 47 Tuc have ax = -3.1 (Grindlay 
et al. 2002). Comparing the corrected source counts in each 
annulus with the predictions of generalized King models with 
ax = -2, -3, and -4, we find a best fit with ax = -2, suggesting 
similarity in source mass to the turnoff mass stars. However, 
the large uncertainties on our source numbers (especially in the 
core) allow consistency with ax = -3 in a sense. We con- 
clude that the small number of sources does not allow a sensi- 
tive test of the radial distribution, and thus of the mass of the 
X-ray sources. Upcoming Chandra observations of Terzan 5 
(PI: R. Wijnands) should identify more X-ray sources, allowing 
a better probe of the luminosity function and radial distribution. 

2.5. Spectral fitting and variability of faint sources 

We extracted the spectra of the detected X-ray sources (ex- 
cept for EXO 1745-248; see section 2.6) from 1" radius circles 
around the wavdetect positions, with background spectra taken 
from 1-2" annuU, using the psextract CIAO script. We utilize 
the recently released ancillary response file correction^ for the 
time- variable low-energy quantum efficiency degradation in all 

^ See http://cxc.harvard.edu/cal/Links/Acis/acis/Cal_prods/qeDeg/index.html. 



spectral fits in this paper. This correction had a very small im- 
pact on our fits (in part because of the large Nh)- We group the 
spectra with at least 10 counts/bin, and fit models in XSPEC 
(Amaud 1996), all using the near-infrared-derived value of pho- 
toelectric absorption (CLG02). The models used were ther- 
mal bremsstrahlung, a powerlaw, and a hydrogen atmosphere 
appropriate for an accreting neutron star in quiescence (Lloyd 
2003). Simple background subtraction is not the preferred way 
to fit spectra from sources with high background, due to de- 
viations from Poisson statistics in the subtracted counts. We 
verified that Sherpa spectral fitting of source and background 
simultaneously gives results within the errors of our XSPEC 
analysis where we are able to use the same models, but note 
that these low signal-to-noise spectra should be taken with cau- 
tion. 

A thermal bremsstrahlung spectrum of kT >,5 keV seems to 
describe most cataclysmic variables in the globular clusters 47 
Tuc and NGC 6397 (Grindlay et al. 2001a,b). We fixed the 
normalization of the hydrogen atmosphere models, set to (r 
[km]/D[10 kpc])^= 132 (appropriate for a 10 km neutron star 
at 8.7 kpc), and the redshift to 0.306. This tests for compatibil- 
ity with a canonical 1 .4 Mq, 10 km neutron star. A hard power 
law component has been found to dominate qLMXB spectra 
above 2-3 keV in some systems, but has not yet been identified 
in globular cluster qLMXB s that have not displayed outbursts 
(see Heinke et al. 2002). We make no attempt to confirm or rule 
out such a hard component in this paper due to the extremely 
high background (from EXO 1745-248) above 2 keV. 

The results of our spectral fitting are shown in Table 2 (er- 
rors are 90% confidence). All of the four sources on the right 
side of the X-ray CMD (W2, W3, W4, W8) are well fit by our 
hydrogen atmosphere model. Alternative models require steep 
powerlaw spectra {a >2) or bremsstrahlung temperatures be- 
low 2 keV, not seen in bright (> lO^'^ ergs s~') CV or millisec- 
ond pulsar systems (Grindlay et al. 2001a, 2001b, 2002, Pooley 
et al. 2002a). Their high luminosities of 3-6 x 10^^ ergs s"' 
and soft spectra rule out flaring coronal events from BY Dra or 
RS CVn binary systems (see e.g. Grindlay et al. 2001a). There- 
fore we conclude that these are probably qLMXB systems. The 
five sources on the left of the CMD (W5, W6, W7, W9, WIO) 
are well fit by bremsstrahlung spectra with kT>5 ke V, or power 
law spectra with a= 1 .0- 1 .6, and are not acceptably fit by hydro- 
gen atmosphere models with the normalization fixed to a 10 km 
radius. We conclude that these are most likely highly X-ray lu- 
minous cataclysmic variables, though MSPs and RS CVn stars 
are (less likely) possibilities. 

We test for variabihty by extracting the events files be- 
tween 0.5 and 2.5 keV within 1" of each source, and using the 
IRAF PROS task vartst. This task applies the Kolmogorov- 
Smirnov and Cramer- Von Mises tests (Daniel 1990) to attempt 
to disprove the hypothesis that the source flux is constant. We 
choose only events below 2.5 keV because the events at higher 
energies are generally dominated by the background from EXO 
1745-248. Only two sources are found to be significantly vari- 
able in both tests, W9 and WIO (identified as probable CVs 
above). Both are variable with 99% confidence. This is con- 
sistent with the trend for qLMXBs of Lx ~ 10^^"^^ ergs s~' to 
show generally constant flux, compared to globular cluster CVs 
that are often highly variable (Grindlay et al. 2001b, Heinke et 
al. 2002, but cf. Becker et al. 2003). 
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2.6. Spectral fitting ofEXO 1745-248 

When observed by the ROSAT All-Sky Survey at a luminos- 
ity of log(Lx)=35.3, the 0.5-2.5 keV spectrum of EXO 1745- 
248 was well fit by a powerlaw of photon index a=1.2, using a 
hydrogen column of Nh = 2 x 10^^ cm"^ (Verbunt et al. 1995). 
A ROSAT HRI observation in March 1991 found it at a lumi- 
nosity of log(Lx)=34.6, assuming the same spectrum (Johnston 
et al. 1995). Only a broad-band X-ray spectrum can provide 
detailed information about the nature of the LMXB, by simulta- 
neously constraining the photoelectric absorption column, any 
iron lines or edges, and the overall spectral shape. Fortuitously, 
there exist RXTE data simultaneous with our Chandra data, 
which we also analyze. 

To study the Chandra spectrum of EXO 1745-248 in out- 
burst, we first selected an annulus from 5 to 9 arcseconds from 
our best-determined position, and removed the portion affected 
by the readout streak. This annulus is not affected by pileup, 
and includes only one faint source (W4) which should not sig- 
nificantly affect the spectral fitting (contributing ~ 17 of 24739 
counts). We extracted the spectrum using bins of width 75 eV 
(oversampling the energy resolution) and using a background 
region on the same chip 2' away. We also extract a lightcurve 
from the same region. We see a gradual rise from 0.4 counts s~' 
to 0.5 counts s"' during the first observation, while the second 
(5 ksec) observation five days later shows a rate of 0.9 counts 
s~' (consistent with the rise in the ASM countrate). Flicker- 
ing is clearly seen, but we put off detailed study of temporal 
variability for future work. An absorbed powerlaw of photon 
index 0.24 provides a good fit to the continuum, with a strong 
broad iron fine (EW AAQ^j^^ eV, d = 0.5 ± 0.1 keV) required 
at 6.7 keV. However, we note that the spectrum will be altered 
both by the effects of dust scattering and by the energy depen- 
dence of the Chandra mirror psf (e.g.. Smith, Edgar & Shafer 
2002). Therefore we extracted the spectrum from the readout 
streak, which does not suffer this hardening since the counts are 
recorded in the image core. 

We used the CIAO task acisreadcorr to identify and repo- 
sition events in the EXO 1745-248 readout streak. We use a 
strip of width 2 pixels and omit a region of 25 pixels radius 
around the source location, to avoid photons from the X-ray 
halo and regions that suffered pileup. Pileup is not an issue 
in the selected region, as the 8750 counts selected give 0.15 
cts (0.841 s CCD frame)"', spread along a 78x2 pixel column. 
We extracted those events using the CIAO script psextract, and 
chose a background region adjacent to the strip (adjusting the 
BACKSCAL parameter by hand) . We ignore data above 1 keV 
and below 0.7 keV (which is almost entirely background pho- 
tons, considering the high absorption to this source). We fit this 
Chandra readout streak spectrum along with the simultaneous 
RXTE data to understand the full spectral shape. 

From the RXTE data, only the PCA STANDARD 2 data 
and the HEXTE Archive data (from both HEXTE clusters 
separately) are analyzed in this paper. We choose time in- 
tervals when three of the PCA units were on (the maxi- 
mum during the observation) and the elevation above the 
earth's limb was greater than 10 degrees. Selecting only the 
top xenon layer data from the PCA, we use PCABACKEST 
version 3.0 (released Feb. 1, 2002) and PCARSP version 
8.0, and correct the PCA data manually for deadtime. We 
used the HEXTE response matrices hexte_97mar20c_pwa.rmf 
and hexte_97mar20c_pwb013.rmf, and corrected the HEXTE 
data for background and deadtime using the HXTBACK and 



HXTDEAD (version 2.0.0) scripts. We eliminate 200 seconds 
around a possible Type 1 X-ray burst, which we do not investi- 
gate in this paper (Type I bursts are common from EXO 1745- 
248; see Inoue et al. 1984). This gave total corrected exposure 
times of 2060, 731, and 777 seconds for the PCA and HEXTE 
clusters A and B respectively. We analyze PCA data from 3 to 
25 keV, and HEXTE data from 25 to 125 keV. Following Bar- 
ret et al. (2001), we add systematic errors of 0.5% to the PCA 
data below 15 keV, and 1% to PCA data above 15 keV using 
the FTOOLS GRPPHA. We leave an overall normalization free 
between the Chandra , PCA and HEXTE data, but link all other 
parameters between them in a joint fit. (The PCA and HEXTE 
relative normalizations are generally not well-calibrated; we 
find the HEXTE normalization to be 40% lower than the PCA 
normalization.) 

The standard models for fitting neutron star LMXB spec- 
tra are an absorbed multicolor blackbody, or a simple black- 
body, with an additional hard component due to comptonization 
of soft photons by hot electrons (generally assuming, as here, 
a spherical geometry); see Barret et al. (2000), and Sidoli et 
al. (2001). Fitting a simpler model to the RXTE data, consist- 
ing of DISKBB plus a gaussian and powerlaw, gave a photon 
index of 1.5, but failed to fit the HEXTE data due to a dropoff 
of flux at high energies (x^=8.1 for 230 degrees of freedom). 
We begin by fitting our RXTE spectra with a model consisting 
of an absorbed multicolor blackbody (DISKBB; Makishima et 
al. 2000), a comptonization model (COMPTT; Titarchuk 1994), 
and the Fe-line gaussian, and with another model consisting 
of a blackbody, COMPTT, and the gaussian. The RXTE fits 
for both models are significantly improved by the addition of a 
smeared iron edge ("smedge"; Ebisawa et al. 1994) near 8 keV, 
as expected for reflection of Comptonized hard X-rays from a 
disk. (An F-test gives a probability of 10"^ that the smeared 
edge is not needed.) For the rest of this analysis, the "standard 
model" shall refer to PHABS(DISKBB + COMPTT + GAUS- 
SIAN)*SMEDGE*CONSTANT The data/model ratio for PCA 
data, fit with the standard model with the normalizations of the 
gaussian and smeared edge set to zero, are shown in Figure 6. 

The Chandra spectrum extracted displays what appear to 
be emission lines at ~1.95 and 2.1 keV. These features are 
not apparent in either the adjacent background, or the annu- 
lus spectrum described above. Similar features are also seen in 
flie readout streak spectrum of GX 13-1-1 (Smith et al. 2002), 
and in the high count rate continuous-clocking spectrum of RX 
J 170930.2-263927 (Jonker et al. 2003). The best-fitting stan- 
dard model to the RXTE data does not give an acceptable fit to 
the Chandra data above 5 keV, where they overlap, nor does 
any other model fit both the Chandra and RXTE data; see Fig- 
ure 7. For the best-fit standard model, the reduced x^=1.41 
for 225 degrees of freedom (null hypothesis prob.=5.7 x 10"^). 
Note that the Chandra spectrum appears to have been shifted 
in energy compared to the model predictions. We understand 
this effect as being due to a gain shift between the caUbrated re- 
sponse of the timed exposure mode, and the actual response to 
events occurring during the readout period, when the voltages 
in the CCD are expected to be different. By measuring the dif- 
ference between the model and data edges near 2 keV, and the 
model (derived from RXTE fits) and data between 5 and 9 keV, 
we estimate the gain shift at 7±1%. 

We attempt to compensate for this gain shift by increasing 
the size of the energy bins in the events file by 7%, to 15.6 
eV from 14.6 eV, while keeping the same response, ancillary 
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response, and background files and extracting the altered-gain 
spectrum in the same way. Fitting this shifted-energy spectrum, 
along with the RXTE spectra, gives a greatly improved fit (see 
Figure 8). The standard model fit to these spectra gives a re- 
duced x^=l.06 for 225 dof (prob=0.250). The derived param- 
eter ranges are robust; spectral fits to RXTE data alone give 
similar parameter values (with larger uncertainties for the soft 
components, and Nh virtually unconstrained). Although this 
simple calibration must be used with a great amount of caution, 
we believe that the quality of the resulting spectral fits supports 
our decision to utilize it for this analysis. 

We list the best-fitting parameters (and 90% confidence er- 
rors) for the two models in Tables 3 and 4, along with the de- 
rived values for the radii of the thermal components and the 
ratio of luminosities in the two components /. We have not 
corrected the multicolor disk-blackbody spectral parameters for 
the spectral hardening factor expected in high accretion rate 
systems (see Shimura & Takahara 1995, Merloni et al. 2000), 
for ease of comparison to the work of SPOOl and Parmar et 
al. (2001) below. The averages of the best-fit luminosities (from 
the two fits' PCA normalizations) are Lx(0.5-10)=2.1 x 10^^ 
and Lx(0.1-100)=6.6 X 10^' ergs s"^ Several authors (e.g. 
Bloser et al. 2000) have commented upon an absolute uncer- 
tainty in the PCA flux normalizations of order 15%, which we 
do not attempt to compensate for. We note that the ranges of A^^ 
required by these fits (1.3-1.9 x 10^^ cm"^) are sUghtly larger 
than the infrared-derived estimate of CLG02, 1.2 x 10^^, indi- 
cating probable internal absorption in the LMXB system. 

3. DISCUSSION 

3.1. Ultra-compact nature ofEXO 1 745-248 

We compare our results from fitting combined Chandra and 
RXTE data with our "standard model" spectrum, to the spec- 
tral fits of SPOOl and Parmar et al. (2001) for the other glob- 
ular cluster LMXBs, observed with BeppoSAX. Five compar- 
isons using the DISKBB + COMPTT model in SPOOl separate 
the "normal" (not ultra-compact) globular cluster LMXBs in 
Terzan 2, NGC 6440, NGC 6441, and Terzan 6 from flie ultra- 
compact (binary periods < 1 hour) LMXBs in NGC 6624, NGC 
1851, NGC 6712, and (probably) NGC 6652. SPOOl note fliat 
these comparisons suggest that the DISKBB model is physi- 
cally meaningful only for the ultracompact LMXBs. 

• The ultracompact LMXBs show an inner DISKBB tem- 
perature below 1 keV, while the normal LMXBs show kTin be- 
tween 1.9 and 3.5 keV. EXO 1745-248 shows kTi„ = 0.80!-^| 
keV. 

• The seed photon temperature kTo for the ultracompact bi- 
naries is roughly equal to the inner DISKBB edge temperature, 
while for the normal LMXBs the seed photon temperature is 
4-5 X lower. EXO 1745-248 seems to have similar values for 
kTo (1.29!;|!,' keV) and kTjn (0.80!;^|* keV). 

• The inferred inner radius Rin (cos i)"^ of the normal 
LMXBs is very small (0.3 to 1 km), while for the ultracom- 
pact LMXBs it is 3 to 45 km. EXO 1745-248's implied Rin 
(cos i)°-5 is 8.7^^:5 km. 

• The seed photon emission radius /?vt' = 3 x 
IQ'^dyJ fcomptt/(^ +y)/ikTa)^ km, where d is the source distance 
in kpc, fcomptt is the luminosity in the Comptonized component 
in ergs cm"^ s"', and y is the COMPTT Comptonization pa- 
rameter, y = AkTeT^ / mgC^. Ry/ is similar to the inner disk radius 
Rin (cos i)°^ in the ultracompact LMXBs, while it is ~ 10-50 x 
larger for the normal LMXBs. For EXO 1745-248 we calculate 



Rw=5A-^J^l km, similar to its Rin (cos i)°^ = 8.7!^;^ km. (We 
also note that y = 5.2^^^ for EXO 1745-248, wifliin the range 
1-7 of most of the LMXBs studied in SPOOl). 

• Finally, the implied mass of the compact object, derived 
from the DISKBB extrapolated bolometric luminosity and the 
inner disk temperature (following Makishima et al. 2000, in- 
cluding spectral hardening), is consistent with 1-3 Mq for the 
ultracompact binaries, but is nearer 0.1 Mq for the normal 
LMXBs. We derive 1.6!jj;^ Mgfor EXO 1745-248. 

Each of these comparisons strongly indicates that EXO 1745- 
248 is an ultracompact LMXB. Some of these spectral differ- 
ences have been used by Parmar et al. (2001) to identify the 
bright LMXB in NGC 6652 as a probable ultracompact sys- 
tem, which is supported by the short optical periods proposed 
by Heinke et al. (2001). The identification of EXO 1745-248 
as a probable ultracompact LMXB makes five of 13 luminous 
LMXBs in globular clusters ultracompact, a far greater pro- 
portion than in the field. This strengthens the conclusion of 
Deutsch et al. (2000) that dynamical effects are probably re- 
sponsible for the generation of ultracompact LMXBs in globu- 
lar clusters. 

Another interesting result is the strength of the ionized iron 
fine in EXO 1745-248, with an equivalent widfli (EW) - 190 
eV. This is many times larger than the average upper limit (25.6 
eV) on the equivalent width of an iron Une for the nine globular 
cluster LMXBs with good statistics reported by SPOOl and Par- 
mar et al. (2001). Only two other globular cluster LMXBs have 
identified iron fines; 4U 1820-30 in NGC 6624 (EW=13^J? eV, 
SPOOl; EW=31!}| eV at 6.6 keV, Asai et al. 2000; EW=27- 
94 eV in numerous RXTE observations, Bloser et al. 2000) 
and Terzan 2 (EW 2\^^ or 27±10 eV depending on contin- 
uum model. Barret et al. 2000). Among these LMXBs, only 
4U 1820-30 in NGC 6624 has a reported reflection edge fea- 
ture (at 7.7 or 8.9 keV; Bloser et al. 2000), in its low or "island" 
state. The strength of the iron line may be increased by the high 
metallicity of Terzan 5, suggested by CLG02 to be greater than 
solar, compared to the significantly subsolar metallicity for the 
other globular s with LMXBs except Liller 1. We note that the 
6.55;';q7 keV energy a\\da = Q3\^W keV breadth of this line are 
similar to those of field neutron star LMXBs reported in Asai et 
al. (2000), 6.56 keV and 0.2±0. 1 keV (0.5 keV FWHM) respec- 
tively, although this line is somewhat stronger than in most field 
LMXBs. The iron Unes studied by Asai et al. (2000) were inter- 
preted as radiative recombination in an accretion disk corona, 
with the breadth being due primarily to Compton scattering, 
with possible small contributions from Doppler shifts (due to 
Keplerian motion near the NS) and line mixing from plasmas in 
different ionization states. The location of our line and smeared 
edge suggest iron that is less highly ionized than He-like Fe 
XXV K lines (6.68 keV line and 8.8 keV edge). Disk reflection 
from carbon-like Fe XXI (6.54 and 8.3 keV line and edge) is 
more consistent with our line and edge energies. Further analy- 
sis of the RXTE outburst spectrum and variability of EXO 1745- 
248 is being performed by Homan et al. (in preparation). 

Kuulkers et al. (2002) have recently presented a BeppoSAX 
spectrum of EXO 1745-248, taken 15 days later during the 
same 2000 outburst. Kuulkers et al. (2002) find an inter- 
steUar absorption of Nh = ^-^-oj x 10^^ cm~^ when using the 
DISKBB + COMPTT model, and Nh = 2.3!g f x 10^^ using the 
BBODYRAD H-COMPTT model. These values are inconsistent 
with our (more precise) values for Nh, suggesting that the ab- 
sorption colunm may have changed during the outburst. Their 
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Nh values are inconsistent with the optically derived extinction 
value (CLG02; see figure 1 in Kuulkers et al. 2002), requir- 
ing an absorption column within the LMXB system ten times 
larger than shown in any other globular cluster LMXB. If this 
were the case, we would expect a high inclination and signifi- 
cant dips due to variable absorption. While these are not seen in 
our Chandra and single-observation RXTE lightcurves, R. Wi- 
jnands notes (priv. comm.) that such dips are indeed seen in 
the full RXTE lightcurves of EXO 1745-248's 2000 outburst 
(J. Homan, in prep.). This may imply that the system is at high 
incUnation, though the observed disk reflection component sug- 
gests a low inclination. We also note that their spectral param- 
eters, while not agreeing with ours in every detail, support our 
claim that EXO 1745-248 is ultracompact, particularly in their 
smaU value of kT,„ for their DISKBB + COMPTT fit. 

3.2. Terzan 5 cluster parameters 

It has long been suspected that Terzan 5 has one of the high- 
est rates of close encounters between stars of any Galactic glob- 
ular cluster. The "collision rate", or rate of close encounters 
given by r oc p^rl ju (where po is the central density, r^ is the 
core radius, and a is the central velocity dispersion), is ex- 
pected to predict the relative rates of formation of accreting 
binary neutron star systems by two-body encounters (Verbunt 
& Hut 1987). Thus, the similar numbers of X-ray sources in 
NGC 6440 and Terzan 5 might seem a surprise, as Terzan 5 
has been predicted to show three times as many collision prod- 
ucts as NGC 6440 and 17% of the total of such objects in the 
Galactic globular system (Verbunt 2002). 

Our calculation of the central density of Terzan 5 uses the 
extinction-corrected central surface brightness /ii/(0) = 20.5, 
combining the star-count profile of CLG02 for the inner core 
with the surface brightness profile of Trager et al. (1995) for 
normalization beyond 10". We use the central concentration pa- 
rameter c=2.0, core radius rc=7"9, heliocentric distance of 8.7 
kpc, and A^/=6.75 from CLG02, as well as My =-7.91 from 
the updated Harris (1996; rev. 1999) catalog. Following the 
prescription of Djorgovski (1993), our result is that Terzan 5's 
central density is 1 .7 x lO^L© pc""'', significantly less than given 
by Djorgovski (1993), Harris (1996), or CLG02. The former 
two studies used a larger value of Ay, which produces a larger 
correction of the surface brightness. CLG02 had scaled the cen- 
tral surface brightness value of Djorgovski (1993) for changes 
in the profile, without updating the value of Ay. We note that 
this calculation is consistent with the lower limit of 5.0 x 10^ 
M0pc~^ derived by Lyne et al. (2000) from the acceleration- 
induced E of pulsar Terzan 5 C, provided that M/L> 3.0. Us- 
ing this central density, and the distance and core radius from 
CLG02, the collision rate in Terzan 5 is similar to that in NGC 
6440 (5.9% of the total collision rate in the globular system), 
instead of three times larger This is consistent with the re- 
sults of our artificial star tests (section 2.3), which suggest that 
Terzan 5 may contain ~1 1 sources in the range 10^^"^^ ergs s~\ 
compared to 1 1 in NGC 6440. A full analysis of the luminosity 
function and density weighting of quiescent LMXBs in globular 
clusters will be presented in Heinke et al. (2003, in prep.) We 
also revise the estimate of the central relaxation time of Terzan 
5 to 2 X 10^ years from CLG02's value of 4 x 10^ years, follow- 
ing Djorgovski's (1993) scaling of trc oc (^Q^r\. This indicates 
that the cluster is not as close to the verge of core collapse as 
suggested by CLG02. The high density of compact binaries in 
the core seems to be due primarily to the high density of the 



cluster's massive core. The large numbers of millisecond pul- 
sars (60-200) in the core of Terzan 5 estimated by Fruchter and 
Goss (2000) would require a long period of high core density 
to form the millisecond pulsar progenitors. MSPs may be most 
efficiently formed early in the cluster history by intermediate- 
mass X-ray binaries (e.g. Davies & Hansen 1998). Therefore it 
seems Ukely that Terzan 5's MSPs were formed early in its his- 
tory, as were those in 47 Tuc (Grindlay et al. 2002), which has a 
roughly similar inferred ratio (~ 1 5 - 50) of MSPs to qLMXBs. 

4. CONCLUSIONS 

We have presented a reasonable infrared candidate to EXO 
1745-248 in Terzan 5, identified by its blue color and posi- 
tional coincidence with the boresighted Chandra position. We 
have assembled a broad X-ray spectrum using a simultaneous 
RXTE observation and the Chandra spectrum from the readout 
streak, shghtly altering the energy scale of the Chandra read- 
out spectrum to account for observed gain variation during the 
readout. We utilized the empirical comparisons of SPOOl to 
indicate that this LMXB appears similar to other ultracompact 
LMXBs in globular clusters, suggesting that EXO 1745-248 is 
the fifth ultracompact LMXB known in a globular cluster. We 
also identify a broad, strong 6.55 keV iron line, the strongest 
(EW=188!83 eV) yet discovered in a globular cluster LMXB, 
with an accompanying smeared ^8.1 keV iron edge. 

The superb resolution of Chandra has allowed us to identify 
nine faint X-ray sources within 30" of an LMXB in outburst 
in Terzan 5. Spectral fitting with bremsstrahlung and power- 
law models, and a neutron star hydrogen atmosphere model 
(Lloyd 2003), suggests that four of these sources are qLMXBs, 
while five are candidate CVs. Artificial point source testing 
suggests that we are missing ~30% of the sources in the range 
Lx(0.5 -2.5 keV) = lO^^"^^ ergs s"^ due to flie presence of flie 
LMXB in outburst. This implies a total cluster population of 
ll-4;l^;8(la) sources with Lx > 10^^ ergs s"' (excluding the 
LMXB). A recalculation of the central density of Terzan 5 from 
updated cluster parameters gives log(/9o)=5.23, suggesting that 
Terzan 5 is not the richest of the globular clusters in stellar en- 
counter products and is not as dynamically unstable as previ- 
ously thought (CLG02). Thus we find that the numbers of X- 
ray sources in Terzan 5 are consistent with the numbers discov- 
ered in other globular clusters and the currently favored forma- 
tion methods. 

Upcoming Chandra observations of Terzan 5 and NGC 6440 
(PI: R. Wijnands) wiU allow us to better constrain the variabiUty 
and spectra of the qLMXBs in those clusters. V and I observa- 
tions of Terzan 5, with the HST Advanced Camera for Surveys 
in HRC mode at times when EXO 1745-248 is in outburst vs. 
quiescence, could unambiguously verify the identification pro- 
posed here. 

C. H. flianks J. M. Miller, R. Wijnands, T. Gaetz, R. K. Smith, 
R. Edgar, and the anonymous referee for useful comments that 
have improved this paper. CH also thanks S. Wachter for ac- 
cess to unpublished data. This work was supported in part by 
Chandra grant GO0-1098A and HST grant GO-7889.01-96A. 
RXTE data and results provided by the ASWRXTE teams at 
MIT and at the RXTE SOP and GOF at NASA's GSFC. The 
Guide Star Catalogue-11 is a joint project of the Space Tele- 
scope Science Institute and the Osservatorio Astronomico di 
Torino. This research has made use of the data and resources 
obtained through the HEASARC on-line service, provided by 
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Table 1 

Names, Positions and Counts of Detected Sources 



Wmiv*/>a ^Janm iX anal I 

oUurLc iconic \i^d.ijd} 


P A 
t\J\ 




u. j-^.j Ke V 


U.J- 1 .J Kc V 


1 < If f»V 
l.J-D Kc V 








(^COUnlS ) 


(^COUnlS ) 


(^COuniS ) 


EXO 1745-248 (LMXB) 


48:05. 196±.015 


46:47.40±.20 








CXOGLB J174806. 1-244642.9 (W2) 


48:06.154±.005 


46:42.68±.07 


37±10 


11±4 


20±10 


CXOGLB J174805.3-244637.7 (W3) 


48:05.370±.004 


46:37.66±.05 


77±13 


30±6 


41±13 


CXOGLB J174804.7-244644.6 (W4) 


48:04.799±.003 


46:44.64±.06 


17±19 


15±7 


-1±21 


CXOGLB J174804.4-244638.2 (W5) 


48:04.402±.006 


46:38. 15±.09 


81±11 


13±4 


80±12 


CXOGLB J174804.3-244703.8 (W6) 


48:04.366±.004 


47:03.75±.07 


180±16 


22±5 


208±17 


CXOGLB J174804.2-244641.8 (W7) 


48:04.214±.005 


46:41.79±.09 


85±11 


10±4 


86±12 


CXOGLB J174804.2-244648.4 (W8) 


48:04.225±.05 


46:48.34±.2 


34±12 


15±5 


30±12 


CXOGLB J174804.0-244640.5 (W9) 


48:04.059±.009 


46:40.53±.06 


50±10 


4±3 


53±11 


CXOGLB J174803.5-244649.2 (WIO) 


48:03.539±.005 


46:49.24±.08 


69±11 


8±4 


75±12 


Sources not associated with Terzan 5 


CXOU J 174803. 3-244854.1 (Star 1) 


48:03. 334±.004 


48:54.04±.06 


45 ±8 


36±7 


9±5 


CXOU J 17475 1.7-244657.4 


47:51.726±.010 


46:57.42±.ll 


36±7 


9±4 


35±7 


CXOU J174814.7-244802.2 


48:14.703±.011 


48:02.21±.12 


29±7 


7±4 


28±6 


CXOU J174812.6-244811.1 (Star 2) 


48:12.649±.010 


48:11. 10±.13 


28±8 


20±6 


6±5 



Note. — Sources detected in and around Terzan 5. RA and Declination values for cluster members begin with 17: and -24:, 
respectively. Errors are centroiding errors within the (shifted to match GSC 2.2) Chandra frame, and do not include systematic 
errors incurred in matching to other frames (absolute errors perhaps 0"5; see text). Counts in each energy band were determined 
by subtracting the average count rate in surrounding 2" aimulus from the counts in 1" circle around source position (leading to one 
negative entry). 



Table 2 

Spectral Fits to Faint Sources 



Source H-atmosphere Bremsstrahlung Powerlaw 0.5-2.5 Lx 

(kT.eV) (x^/dof) (kT.keV) (x^/dof) (a) (x^/dof) (ergs s'^) 



W2 


98!i 


1.19/14 


0.7!°:^ 


1.16/13 


3 Q+l-l 


1.31/13 


4.1x1032 


W3 


104+5 


1.36/22 




1.05/21 


1 1+0.8 
■'■^-0.6 


1.20/21 


5.6x10^2 


W4 


90!^, 


1.37/11 


5+2-5 
^•-'-0.5 


1.49/10 


4 

^•"-2.6 


1.5/10 


2.7x10^2 


W5 


92 


3.0/16 




1.29/15 


1 7+0.5 
^•'-0.5 


1.34/15 


1.9x10^2 


W6 


100 


4.42/33 


>29 


1.08/32 




1.05/32 


3.5x10^2 


W7 


90 


2.50/16 


>6 


0.51/15 


1 O+0.5 
^ -^-0.5 


0.51/15 


1.5x10^2 


W8 


92!^ 


0.84/19 




0.78/18 


1 a+1.8 


0.82/18 


3.1x10^2 


W9 


79 


2.27/13 


>2 


1.37/12 


1 O+0.9 
^■•'-0.8 


1.37/12 


8.6x10^1 


WIO 


88 


2.92/15 


>7 


0.95/14 


1 1+0.6 
^•^-0.6 


0.94/14 


1.3x10^2 



Note. — Spectral fits to faint cluster sources, with background subtraction, in XSPEC. All fits 
include photoelectric absorption fixed at the cluster Nh, 1.2x 1022 cm"2. Hydrogen atmosphere fits 
are made with radius fixed to 10 km. X-ray luminosities are unabsorbed for the range 0.5 to 2.5 keV, 
from hydrogen-atmosphere NS fits (W2, W3, W4, W8) or thermal bremsstrahlung fits. Errors in all 
the tables are 90% confidence for a single parameter. 
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Table 3 

Spectral Fits to LMXB EXO 1745-248 Continuum 



Model 


Nh 
(10^2 cm-2) 


(kTo) 


COMPTT 




BB or DiskBB' 
(kT) (R, km) 


(f*) 


Xl/dof 


DBB 
BB 


1.72+'20 


1.29!f5 

1 21 +-24 
^•^^-.09 


10.2!-5 
10.2^5 


8.i!i 
B.i^i 


0.80+?^ 




0.12 
0.07 


1.06/225 
1.06/225 



Note. — Spectral fits to the LMXB EXO 1745-248 from RXTE data and the Chandra readout streak 
with energy scale alteration (see text), with background subtraction in XSPEC. A gaussian, a smeared edge, 
photoelectric absorption, and comptonization using the XSPEC model COMPTT is included in both fits. 
The first also includes a multicolor disk-blackbody (DBB), while the second includes a blackbody (BB). 
Values are in keV unless otherwise noted (f and Tp are unitless). Values from normalizations utilize the PC A 
normalizations; uncertainty in PCA absolute normalization is not included. 

^kT and R refer to the inner edge of the diskbb for model 1 (where R is R,„cos(i)''-^), and to the blackbody 
for model 2. 

''f is the ratio of the blackbody or diskbb flux to COMPTT flux, over 0. 1- 100 keV range. 



Table 4 

Spectral Fits to LMXB EXO 
1745-248 Features 

Gaussian 



kT 6.55!-^^ 
EW (eV) 188!^^ 

Smeared edge 
IT 8.1±0.8 

'max yJ-^-J-i^ 

Width 2.6!^:^ 



Note. — Spectral fits to the 
Fe-line gaussian and disk reflection 
smeared edge in the spectrum of 
the LMXB EXO 1745-248. The 
spectral fits in this table are to the 
"standard model" including a multi- 
color disk-blackbody (see text); line 
and edge parameter values for the 
other model considered in the text 
are within the errors in this table. 
Values are in keV unless otherwise 
stated, except Tmax which is unitless. 
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Fig. 2. — Chandra ACIS-I image of the globular cluster Terzan 5, energy range 0.5-2.0 keV. The dominant feature is the pUed-up halo from the LMXB EXO 
1745-248 in its high state. The streak is due to out-of-time LMXB events recorded during the frame transfer, and the position of the LMXB has no good events due 
to pulse saturation. Several other sources are visible, marked with 1" circles and indicated with their shorthand names. The cluster core is indicated with a 7.9" (1 
rc) circle. 




Fig. 3.— Finding chart, showing stars A and B in FllOW (upper left) and F187W (lower left) //5r NICMOS data, and the ratio image F110W/F187W (upper 
right). The best estimate of the uncertainty in the position of EXO 1745-248 is indicated by the 0.8" (2 ct) error circle. Note that star A is the bluest object in this 
field (some red giant cores appear blue due to saturation). The results of our manual psf fitting of stars A and B are shown at lower right. 
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12 3 

F110W-F187W 



Fig. 4. — Aperture photometry color-magnitude diagram of Terzan 5, from the HST NIC2 camera (using highly dithered imaging via the drizzle algorithm) on the 
core of Terzan 5. Stars from the photometry of EGCOl within 1" of the seven X-ray error circles are circled. Xs mark the locations of variable stars from EGCOl 
(X2 is the dimmer one). Two stars that are more than 3 cr to the blue of the distribution are indicated by < . >. The lower is the blend of stars A and B (see text). 
The separately derived positions of stars A and B are indicated, with errors from the photometry. The horizontal branch and the expected blue straggler sequence, 
terminating in a Kurucz model (Kurucz 1992) of a 1.6 Mq star (twice the tumoff mass) are indicated. 



16 



Heinke et al. 




-10 12 

Xcolor = 2.5 log([0.5-1.5keV]/[L5-6keV]) 




-10 12 

Xcolor = 2.5 log([0.5~1.5keV]/[L5-6keV]) 



Fig. 5. — X-ray color-magnitude diagrams of Terzan 5. Upper panel produced from the absorption-corrected counts in the 0.5-1.5, 0.5-4.5, and 1.5-6 keV energy 
bands (lower and left axes; upper and right axes provide the observed colors and magnitudes). We correct roughly for photoelectric absorption by shifting the data 
-hO.57 units on the left axis and -1-1.86 units on the bottom axis. However, note that the effects of absorption are not uniform. The probable qLMXBs are marked by 
crosses, and probable CVs by filled triangles. We plot representative errors for a few points. Background subtraction leaves W4 with -1 counts in the hard band; we 
show its Xcolor lower limit and counts range (it is only clearly detected in the soft band). The Xcolor of EXO 1745-248's X-ray halo is indicated with a dotted line 
and labeled 'LMXB Xcolor' . Lower panel produced using the same absorption-corrected color for x axis, but using the vinabsorbed 0.5-6 keV luminosity for y axis. 
Luminosities derived from best-fitting hydrogen atmosphere fits for W2, W3, W4, and W8, and best-fitting thermal bremsstrahlung fit for W5, W6, W7, W9, and 
WIO. Luminosity errors are based on 0.5-4.5 keV counting statistics, and do not include imcertainties in spectral fitting. 
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Fig. 6. — Ratio of EXO 1745-248 PCA data to standard model (see text), with normalizations of the gaussian iron line and smeared edge set to zero to show the 
relative contributions. 
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Fig. 7. — EXO 1745-248 spectrum from Chandra readout streak data, RXTE PCA, and RXTE HEXTE data, simultaneously fitted with the standard model (see 
text). Note the poor fit to the Chandra data near 2 keV and the pronounced wave above 5 keV, suggesting an incorrect energy calibration. 
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Fig. 8.— EXO 1745-248 observed (top) and unfolded (bottom) spectra, with best-fit absorbed (DlSKBB-l-COMPTT-l-GAUSSlAN)xSMEDGE model spectrum. 
Here the energies assigned to Chandra readout streak events have been increased by 7%. Relative normalizations of the different instruments are allowed to vary 
(see text). The 6.55 keV iron line is visible in both the RXTE and Chandra unfolded spectra, along with the DISKBB and COMPTT components (the COMPTT is 
the harder component). 



